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Abstract 

The  singular  stress  problem  of  a  peripheral  edge 
crack  around  a  cavity  of  spherical  portion  in  an 
infinite  elastic  medium  when  the  crack  is  sub¬ 
jected  to  a  known  pressure  is  investigated.  The 
problem  is  solved  by  using  integral  transforms 
and  is  reduced  to  the  solution  of  a  singular  in¬ 
tegral  equation  of  the  first  kind.  The  solution 
of  this  equation  is  obtained  numerically  by  the 
method  due  to  Erdogan,  Gupta  ,  and  Cook,  and 
the  stress  intensity  factors  are  displayed  graphi¬ 
cally. 

Also  investigated  in  this  paper  is  the  penny¬ 
shaped  crack  situated  symmetrically  on  the  cen¬ 
tral  plane  of  a  convex  lens  shaped  elastic  mate¬ 
rial. 

Key  words:  cavity  of  spherical  portion/  periph¬ 
eral  edge  crack/penny-shaped  crack  /SIF. 

1. Introduction. 

The  problem  of  determining  the  distribution 
of  stress  in  an  elastic  medium  containing  a  cir¬ 
cumferential  edge  crack  has  been  investigated  by 
several  researchers  including  the  present  author. 
Among  these  investigations,  the  notable  ones 
are  Keer  et  al.[l,2],Atsumi  and  Shindo[3,4],  and 
Lee [5, 6].  Keer  et  al. [1]  considered  a  circumfer¬ 
ential  edge  crack  in  an  extended  elastic  medium 
with  a  cylindrical  cavity  the  analysis  of  which 
provides  immediate  application  to  the  study  of 
cracking  of  pipes  and  nozzles  if  the  crack  is  small. 


Another  important  problem  involving  a  cir¬ 
cumferential  edge  crack  is  that  concerned  with  a 
spherical  cavity.  Atsumi  and  Shindo[4]  investi¬ 
gated  the  singular  stress  problem  of  a  peripheral 
edge  crack  around  a  spherical  cavity  under  uni¬ 
axial  tension  field.  In  more  recent  years,  Wan 
et  al.  [7]  obtained  the  solution  for  cracks  emanat¬ 
ing  from  surface  senri-spherical  cavity  in  finite 
body  using  energy  release  rate  theory.  In  previ¬ 
ous  studies  concerning  the  spherical  cavity  with 
the  circumferential  edge  crack,  the  cavity  was  a 
full  spherical  shape.  In  this  present  analysis,  we 
are  concerned  with  a  cavity  of  a  spherical  por¬ 
tion,  rather  than  a  full  spherical  cavity.  More 
briefly  describing  it,  the  cavity  looks  like  a  con¬ 
vex  lens.  Here,  we  employ  the  known  methods  of 
previous  investigators  to  derive  a  singular  inte¬ 
gral  equation  of  the  first  kind  which  was  solved 
numerically,  and  obtained  the  s.i.f.  for  various 
spherical  portions.  It  is  also  shown  that  when 
this  spherical  portion  becomes  a  full  sphere,  the 
present  solution  completely  agrees  with  the  al¬ 
ready  known  solution. 

2. Formulation  of  problem  and  reduction  to 
singular  integral  equation.  0 

We  employ  cylindrical  coordinates  (r,  0,  z) 
with  the  plane  2  =  0  coinciding  the  plane  of  pe¬ 
ripheral  edge  crack.  The  spherical  coordinates 
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(p,  9,  (f>)  are  connected  with  the  cylindrical  coor¬ 
dinates  by 

z  =  pcos6 ,  r  =  psin9. 

Spherical  coordinates  (£,  d,  4>)  whose  origin  is 
located  at  z  =  —  5,  r  =  0,  and  is  the  center  of 
the  upper  spherical  surface,  are  also  used.  The 
cavity  is  symmetrical  with  respect  to  the  plane 
z  =  0. 

The  crack  occupies  the  region  z  =  0,  1  < 
r  <  7.  So  the  radius  of  the  spherical  cavity  is 
Co  =  v+++. 

The  boundary  conditions  are: 

On  the  plane  2  =  0,  we  want  the  continuity  of 
the  shear  stress,  and  the  normal  displacement: 


uz(r,0+)  -uz(r,0  )  =  0,  7  <  r  <  00,  (2.1) 

o’rzi.i",  0+)  —  arz(r,  0“)  =  0,  1  <  r  <  00. 

(2-2) 

And  the  crack  is  subjected  to  a  known  pressure 
P(r),  he., 

<fzz{r,  0+)  =  -p(r),  1  <  r  <  7.  (2.3) 

On  the  surface  of  the  spherical  cavity,  stresses 
are  zero: 


crCc((o+)  =  0, 

(2.4) 

cc4Co+)  =  0. 

(2.5) 

We  can  make  use  of  the  axially  symmetric  so¬ 
lution  of  the  equations  of  elastic  equilibrium  due 
to  Green  and  Zerna  [8]  which  states  that  if  ip(r,  z) 
and  ip(r,  z)  are  axisymmetric  solutions  of  Laplace 
equation,  then  the  equations 


2  p'Uj'f' 


dp  di/j 

dr  dr  ’ 


(2.6) 


2puz  =  ^  +  z^  -  (3  -  4u+,  (2.7) 

where  p  is  the  modulus  of  rigidity  and  v  is  Pois¬ 
son’s  ratio,  provide  a  possible  displacement  field. 
The  needed  components  of  stress  tensor  are  given 
by  the  equations 


O’rz  — 


d2cj) 

drdz 

d2(f> 


d V  _  n  _  9  sd± 
Z drdz  ^  ^  dr ’ 


_i_  d2^  on  +Y> 

a~=d^  +  z^-2{1-u)Tz- 


(2.8) 


(2.9) 


The  functions  <+)  and  <+)  for  the  regions  z  >  0 
and  z  <  0,  respectively,  are  chosen  as  follows: 


poo 

<++’,  z)  =  {2u  -  1)  /  ^A^)  J0(+)e^2d£ 

Jo 


+£■ 

n— 0 


Pn( COS  6) 
ln  P-+1 


(2.10) 


poo 

</>( 2){r ,  z)  =  (2v  -  1)  /  £_1A(£)  J0(£r)e^d£ 


tPn(cos/9) 


nn-\- 1 


(2.11) 


71=0 


Here  the  superscripts  (1)  and  (2)  are  taken  for 
the  region  z  >  0  and  z  <  0,  respectively.  The 
functions  +1)  and  +2)  are  chosen  as  follows: 

coo 

ip^(r,z)  =  A^)J0(^r)e~^zd^ 


Pn(cos9) 


n= 0 


P 


,71+1 


(2.12) 


poo 

ip(2\r,z)  =  —  /  A(£)J0(fr)e?2d£ 


+YM-V 


Pn  (cos  6) 


71—0 


P 


,71+1 


(2.13) 


Then  we  can  immediately  satisfy  condition  (2.2) 
by  this  choice  of  functions  (2.10)-(2.13). 

Now  the  condition  (2.1)  requires 

poo 

/  +£)  J0(£r)(f£  =  0,  r  >  7.  (2.14) 


Equation  (2.14)  is  automatically  satisfied  by 
setting 


A(0  =  /  tg(t)Ji(£t)dt. 


(2.15) 
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Then  from  the  boundary  condition  (2.3),  we 
obtain 


roo  00 

/  ~  Va2n 

Jo  n=o 


(2n  +  l)2i32n(0) 
^»2n+3 


1  1  dcj)  d2^ 

~  TA  TTF  +  COS  V- 


(  dQdd  C2  dd 


d(dd 


on  \Vi,  ^2n+l(0)  (  ^ 

-2(1  -  v)  2^  Wi  r2n+3  =  -p(r), 


+(1_2„)sin,|_2(1_,)^g.  (3.2) 


n=0 


To  satisfy  boundary  conditions  on  the  spherical 
surface,  it  is  needed  to  represent  </),' 0  in  (2.10)- 
(2.13)  in  terms  of  £,  1?  variables.  To  do  so  we  uti- 
1  <  r  <  7,  (2.16)  }jze  ^he  following  formula  whose  validity  is  shown 

where  prime  indicates  the  differentiation  with  re-  in  the  Appendix  1.  An  expression  useful  for  the 


spect  to  the  argument. 

By  substituting  (2.15)  into  (2.16),  it  reduces 


to 


present  analysis  is  the  following 
Pn  (cos  6) 


7 r 


tg(t)R(r,  t)dt  -  V  r  (2n+3){T2n(0) 


■)"+ 1 


E 

k= 0 


n  +  k  \  Pn+k{ cost?)  fc 

k  )  (n+k+ 1  ' 


n=0 


(3.3) 


x(2n+  l)2a2n  +  a&2n+l-P2n+l(°)}  =  -p(r); 

1  <  r  <  7,  (2.17) 

where 


Pn(cos9) 


R{r ,  i)  = 


1 


,2  _  ^2 


£|  - 


1 


t  r2  —  t2 


t  >  r, 


E  -  - -A  - 


Thus 

OO 

X] Un  0n+l 

n=0  P 

Pn+fc(cosi?)  , 


Xa-X 

n=0  fc=0 


n  +  k 
k 


^■n+fc+l 


sk  =  ^2  Pj(cos 


rt  \r 


where 


r  >  t. 
(2.18) 


j 


j= 0 


j! 


(■f+i 


3  > 


K  and  E  in  (2.18)  are  complete  elliptic  integrals  Also 
of  the  first  and  the  second  kind,  respectively  and 
a  =  2(1  -  u). 

The  solution  will  be  complete,  if  the  condi¬ 
tions  on  the  surface  of  the  spherical  cavity  are 
satisfied. 


A  .  —  y _ 

J  ^(j-n)!n!' 


b3~ 

iV 


(3.4) 


(3.5) 


roc  /»7 

J  ClA(QMtr)e~^Zd£  =  J  tg(t ) 

r  00 

X  /  (3. 

Jo 


If  we  make  use  of  the  formula  in  Whittaker  and 
Watson[9,pp. 395-396] 


ex' 


p{— £(z+ix cos u+iy sin u)}du  =  2ire  ^2</o(£0> 


3. Conditions  on  the  surface  of  the  spheri¬ 
cal  cavity. 

Equation  (2.17)  gives  one  relation  connecting 
unknown  coefficients  an  and  bn.  The  stress  com¬ 
ponents  besides  (2.8)  and  (2.9)  which  are  needed  aS)  ^  we  are  usinS  the  shortened  notation 

for  the  present  analysis  are  given  by  the  following 
equations 


to  the  inner  integral  of  (3.6),  it  can  be  written 


(5  =  z  +  ix  cos  u  +  iy  sin  u, 


d2<i>  ,  ,  on  x  .d'l/j 

a(C  =  q^2  +  ^cos^qa2  -2(1-u)cosi?— 


then 


+2z/ 


sin  d  dip 

~^~dd' 


(3.1) 


Jo 

l 


C1M£r)Ji(&)e-*zdZ 

1  Ji  (£t)  e~^d£du 


r7r  roo 


—7 r  J  0 
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2?r  J_n  (3  + 


du 


1 


1 

20 


2vr  . 

P 

t 


P 


=  -7T-  /  T~  \  1  +  LT  \du 


P2 


t 2 


g(4)n(-i)V/3x2n 


n=0 


n\ 


du. 


k= 0 


E  ( 2k )  (-^"W 


=  2ir(n Pn(cos  i?) . 

If  we  utilize  (3.8)  and  (3.9)  in  (3.7),  it  can  be 
written  as 

POO 

/  r1«/o(er)J1(et)e-^de 

Jo 

=  g(4)n(- 1)"  1 

72— 0 


n! 


E  (  t"  )  (-i)2”-‘C‘P„(costf) 


+5  —  ^  Pi  (cos  I?). 


=  Ep^cos^) 


k=0 


k  i  P 


/7  77 

tg{t)dt  —  £Pi(cos  i?)  J  g(t)dt , 


where 


$fe  =  -(a-l)  E 


(-l)n(-ir 


[  0-^  <3'+ 

and  [(k+l)/2]  is  the  greatest  integer  <  (fcH- 1) /2. 

It  is  also  necessary  to  express  in  terms  of 
spherical  coordinates  (£,$).  Now,  as  in  (3.4) 


(3-7) 

As  \/3 1  <  |f|,  the  above  series  expansion  is  valid. 
Next, 

/ 32n  =  (-6+Z+ix  cosu+iy  smu)2n  =  (~5+P')2n 

2n  /  § »  x 

\2n—kfa/\k  (3  8) 


Pn  (COS  6) 

nn+ 1  —  2^ 


77.— 0 

where 


l=o 


Pj  (cos  t?) 
£7+1 


Pi,  (3.12) 


p,  =  E  7- 

(j  —  n)!n! 


n=0 


where  (r,  0,  Z)  is  the  cylindrical  coordinates  sys¬ 
tem  centered  at  (r,z)  =  (0,-5).  We  have 
the  following  formula  from  Whittaker  and  Wat¬ 
son^,  p.392] 

P7T  P7T 

/  ( /3')ndu  =  /  (Z  +  ixcosu  +  iysmu)ndu 


We  first  express  following  integral  in  (p,  8)  co¬ 
ordinates 

/»oo  p"y 

J  0P)Jo^r)e~izdi  =  J  tg(t)dt 

POO 

x  /  (3.13) 

Jo 

The  inner  integral  on  the  right-hand  side  of 
(3.9)  (3-13)  is 


1  poo  a  ] 

M0)e~^-Jomdi=- 


t 


dZ 


t 


{rJi(£r)  +  zJ0(£r)}e  ?2J0(£i)d£. 

(3.14) 


Using  the  equation  in  Erdelyi  et  al.[10] 


Jo(&)  =  ~ 


sin  (£x)dx 


(3.10) 


Thus  finally,  from  (3.4),  (3.6)  and  (3.10),  (j)^ 
can  be  written  in  terms  of  spherical  coordinates 
(C,  $)  as 


^  Jt  Vx2  —  t2  1 

equation  (3.14)  is  equal  to 

i  12  r°°  dx 
t  +  t7Tjt  y/X2  -  t 2 

POO 

x9=  /  {r  Ji(£r)  +  zJo(£r)}e~^z+lx^d(, 

Jo 


-i +ii  r  1  - 


—  IX 


t  t  n  Jt  yjx2  —  t2  0r2  +  (z  +  ix )2 


dx 


1  1  2  r°°  1 

t  t  7T 


dx 


1  1 
t  t 


t  Vx2-t2  -2£*cos0  +  (fi)2 

(*)r 


oo  /  N  2n+l 


n=[(fc+l)/2] 


n\ 


E 

71—0 


(-l)nP2ri+1(cos0) 


^2>n 
n\  ’ 

(3.15) 


°@IJTSRD —  Available  Online@www.ijtsrd.com — Volume-2 — Issue-3 — Mar-Apr  2018  Page:985 


INTERNATIONAL  JOURNAL  OF  TREND  IN  SCIENTIFIC  RESEARCH  AND  DEVELOPMENT(IJTSRD)ISSN:2456-6470 


where  we  used  the  generating  function 

1  ~ 


-A 


k-\- 1 ' 


( k  +  2  )(/c  +  3) 


y/l  —  2xcos  6  +  x2 


Pn(cos  6)xn. 


c, 


fc+4 


n=0 


+ 


To  express  (3.15)  in  the  (£,  $)  spherical  coordi¬ 
nates,  we  use  the  following  formula  whose  valid¬ 
ity  is  shown  in  the  Appendix  2. 

2n+l 


Bk  (A  +  1 ) { 2  —  a  —  (A  +  1)(A  +  4)} 
2k  +  1 

-Bfc+2  (A:  +  2)(A  +  3)  (2a  +  A:  +  2) 


/~k-\-  2 
SO 


p2rl+1P2n+l(cOsd)  =  ^ 
fc=0 

\2n+l— k/-k 


(2n+  1)! 
A!(2n  +  1  —  A;)! 


/-fc+4 

so 

1/ 


2A  +  5 

-A;(A;+l)Cofc-^fe+i-Co"i^fc 


k — i it.  A; (A;  +  1)(A  +  1  —  2a) 


x(Syn+1-«(kPk(cos'&).  (3.16)  -Co 

0  If  we  use  (3.16)  in  (3.15),  then  with  (3.12),  we 
get  in  spherical  coordinates  (C,  i?) 


k-\-l 


fc+2~ 


2A  +  1 

(A  +  2  ){(A  +  2)  (A:  —  1)  +  a  —  2} 


=  0. 


V;(1)  =  5ZPfe(cos?9)(^r+^fcl^A:)  +/ 

(3.17)  =  - 

(-i)nG) 


2A  +  5 

(3.196) 

Thus  if  we  multiply  (3.18)  by  k  +  2  and  subtract 
(3.19b)  from  the  resulting  equation  we  find 


Co+2(2  A;  +  1) 


where 


a(2A  +  3)  +  2A(A  +  1) 


fc(2fc+3)C0fc“1<hfc+i 


**  =  -  E 

n=[k/ 2] 

x(-<5) 


(2n  +  1)! 


n! 


2n+l— /c 


A!(2n  +  1  —  A’)! 

r7 


J i  f2^1 

Then  if  we  substitute  these  values  of  and 
given  by  (3.11)  and  (3.17)  into  (3.2),  and 
simplify  the  results  by  using  the  properties  of 
Legendre  polynomials,  from  the  condition  that 
on  the  spherical  surface  C  =  Co>  the  shear  stress 
a(d  =  0,  we  obtain  the  following  equation 

Bk  a  —  (k  +  l)2 


-A 


fc+i 


A  +  3 

•*-&+ 4 


+ 


f + 2  2k  +  1 


j-k—i  t  (2A  T  3)A(A  +  1  —  2a)  T  . 

+Co  k~ - 2A+1 - L+Co+  *k+2k(k+2) 

(3.20) 

If  we  solve  (3.20)  for  6j  using  the  theorem 
which  is  in  the  Appendix  3  and  the  relation, 

,  a  ~  1  T 
$t+1  -  Hi*1’ 

we  find 


/-K- 1-4  /-K 

so  so 

Bk+ 2  (A  +  3)  (2a  +  A  +  2) 

/-k+ 4 

Ao 


+ACo  1<hfc+i  +  Co 


2A  +  5 

A(A  +  1  —  2a) 


where 
IVi(A)  =  - 


(3.21) 
—a  +  A2' 


2A  +  1 


T  A  +  2  +  a  —  (A  +  2)  (A  +  3) 

- 2FT5  —  = 

(3.18) 

Likewise,  from  the  condition  a^( Co,i?)  =  0, 
we  get  following  two  equations, 


,  2  Bi2(2a  +  1)  1T  ,  . 

£3  3  ^  d7 1  (ck  4)  —  0, 


Cq  +i(2A  +  3)A 

a(2A  +  3) +  2A(A  + 1)  V  A+l 

C02fe+3(2A  +  l)A(A  +  2) 


a(2A  +  3)  +  2A(A  +  1) 
Equation  (3.21)  can  be  written  as 


_  ‘  i!(-l)‘-*5* 

0i  ~  2^  u 


9(t ) 


k=0 


(i  —  A)! A!  5k  J 1  t 


(3.19a) 


x  Et^WAOw 


77-— 0 


1  <5  .  , 


(3.22) 
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where 


fk(c,x)  = 


Thus 


(-5) 


i-fe 


£ 


(-1)V)<  (2<+l)!i2>  ”  (2i  +  l)2P2i(0) 


with  (c)g  =  c(c  +  l)(c  +  2)  •  •  •  (c  +  £  —  1). 

Now  from  (3.19a,b)  we  have 

Ak+  i  =  Mi(k)^>k  +  M2{k)'S>k+ 2  + 


where  0 


Mi(ife)  =  - 


A)  =  il^l  +i2^3, 

Cok+3H-a  +  k 2) 


(fc  +  2)(fc  +  3)(2A:  +  1) 

(2fc  +  3)5(/c) 


h(k) 


+  1 


M2(k)  = 


Co2fc+5[  $(*),„  ,  N  .  £l(k  —  2) 


/c  +  3 


fc(/c+l)+: 


2/c  +  5 


(A:  +  2)(2fe  +  7){2a  +  k  +  2){-a  +  (k  +  2)2}-1 
(2k  +  5)h(k  +  2) 

Clk+\k  +  2  )(fe  +  4)  (2a  +  fe  +  2) 


+ 


M3(fc)  = 


A,  =  a-4- 


A2  = 


h(k  +  2) 

5(2a  +  1)(— a  +  1) 

W) 

3(2a  +  l)Cp5 
h(l)  ’ 


with 


0,2  i  =  (+^1  +  A2i&z)5 


■2i 


2i— 1 


+  £ 

k= 0 


(2i)!(-tf) 


2j— fc— 1 


(2*-l-jfe)!(ife  +  l)! 


m, 


where 


m  =  - 


g{t) 


Mi(A:)/fc(  ^  )  +M2(fc) 


1  5 


2  t 


X  fk+2  (  77,  T  +  M3(k)fk+ 4  7 


1  5 


2  t 


dtt. 


J2a2i~ 

i= 0 


r*22+3 


ao 

^3 

\2 


+  J^(ii4'l  +  I24'3) 


2=1 

Wl' 


x<5 


2i(2z  +  l)2(-l)Hi)i 


|— 3  ^  f 


+  g(2z+l)2(-l)++ 


2=1 


22-1 


k= 0 


(2i  —  1  —  k)\(k  +  1)! 


If  we  interchange  the  order  of  summation  in  the 
last  term  of  the  above  equation,  and  use 


'f'  3  =  — 


n 


9(t) 


1 1  t 

r  9(t) 


' 1 


t 


/‘i  bbdt 

ai  ihdt' 


we  finally  obtain 


X+2i 

i= 0 

where 


(2i  +  1)2P2,.(0) 


r»22+3 


=  -/  tg(t)Ti(r,  t)dt, 


Ti(r,t)  = 


r^t2 


OO  2 


££  Mm+\{k')  fk+2m  (  2  ’  t  )  (  -r 


Lfc=0  m=0 


1  5 


1  5 


+  M1/1U^)+^2/3(^? 

with 


l+j(  - 


£+)  =  2  —  a  —  (fc  +  1)(A:  +  4), 
h(k)  =  a(2k  +  3)  +  2k(k  +  1). 

In  a  similar  way  a(s  can  be  found  from  these 
equations  as  follows: 


2’ t 

(k  +  l)\ 

g  {2i  +  l)\-m\U2i)\(5-2i 


h,\  -  = 


i=[k/ 2+1] 

fS\  00 

T 

Also, 


(2*  -  1  -  k)\i\ 


v(2i  +  l)2(-l)^)i^iW 

n\ 


2=1 


v,  ^+i(0) 
Z^2^1  r2i+3 

i=0 


2^ 


i=0 


^1^22+3 


22+1 

£ 

k=0 


2i-k-\-l 


(2i  +  l)l(-6) 

(2i  +  1  -  Jfe)!A:! 
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g(t) 


t 


Nn+l(k)fk+2r 


n= 0 


1  5_ 
2’ t 


dt. 


1 


OO  1 


T2(r’  *)  =  ^3  Nn+l(k)fk  (  - 


fc=0  n=0 


3  5 


2  r 


X  fk+2r, 


1  5 

2’ t 


~  tg(t){R(r,t)  +  S(r,t)}dt  =  p(r), 

7 r  ./i 


1  <  r  <  7, 


where 


7T 


2k  ' 


|  = 

oo  2 


1  A  1 


(fc  —  1)!  t2k  2  ’ 


1  5 


-4  =  Jim  H  Mm+l{k)fk+2m  (^^)hk(- 


k= 0  m=0 
oo 


2  t 


^(-l)k(lH2k  +  lf 


k= 1 


k\r2k 


1  r2fc-lf(4fc+l)G(A:) 
(2k  +  l)(2k  +  2)  |_4fc  —  1  \  h 

l/’  l 


1  j  (fc-1)!  t2fc-2 

+(2/c  +  1)(— §tS-(2 k  -  1)2 k 


H(k) 

(2k  +  l)(4fc  +  5)(2a  +  2fc  +  1) 
H(k  +  l)(4k  +  3) 


F(fe  +  1) 


4A;  +  3 


G(A;-l)\(-l)fc(^  1 


k\ 


t2k 


If  we  change  the  order  of  summation  in  the  above 
equation,  we  get 

!>«%#  = -/Vm  <r,w, 
i=0  Jl 

where 


+ 


(2k  +  l)2(2fc  +  2)(2  k  +  3)  (2a  +  2k  +  1) 
H(k+  1) 


(A’  +  l)!  t2k+2 


where  0 


H(k)  =  a(4k  +  1)  +  4k(2k  -  1) , 
G(k)  =  2-a-2k(2k  +  3), 
F(k )  =  —a  +  (2k  —  l)2. 


OO  1 


B  =  J™  E  E  "»+•  «  A  (  V 7 )  A+2»  ( v  7 


Thus  finally  the  singular  integral  equation  (2.17) 
becomes 


k= 0  n= 0 
oo  /  iU/3\ 


3  (J 


1  <5 


E 

k= 1 


(-!)*(§)* 


2  r 

(2fc+  l)(4A;  +  5) 
2(fc  +  l) 


2  t 


F(k  +  1)  (-l)fc(|)fc  1  (4fc  +  3)(2/c  +  l)(2fe  +  3) 


(3.23) 


H(k  + 1)  k\  t2k 

(-l)k+\\)k+1  l 

(k  +  1)!  t2k+ 2 


//(A;  +  1) 


1 1  1  / 

c;  Fh  j-L  ^  l 

J  2H(1)\ 

01  V  J-J  i  ^  | 

5(r,t)  =  ~-{Ti(r,t)  +  <xT2(r,f)}. 

When  6  =  0,  we  briefly  show  that  this  equation 
completely  agrees  with  what  Atsurni  and  Shindo 
obtained.  Using 

fS\  (~l)k(\)k(2k  +  l)2  1 

- 


Thus 


A+aB  =  t^)k  (!>*(“  +  » 


k= 1 


k\r2k 

k(k  +  l)F(k) 


(—l)k(2k  —  l)(2k  +  1)  f 
H(k)k(k+1)  \ 

'  ‘  {(2k-l)(4:k  +  3)2kG(k) 


+(-l) 


t2  4(4 k  +  3) 
-H(k)G(k-  1)}}  (a 
k  +  1)(2A:  +  1) 


(h)k- i 


(A:  —  1)!  f2fc  2 


/7(A;+  1)(2A’  + 2) 


(2  A:  +  3)(2A;  +  2)=j 


+F(k+1)- 


4:k  +  5  1  (^)fc  1 


4A:  +  3  j  A:!  t2k 


+ 


a  (9 
2/7(1)  \  f2  ~ 


5F(1) 


E 

k= 1 


2(2k  +  l)2 
//(fc)(fc  +  l) 


fc(/t  +  l)F(A-)  - 


f2  4(4 k  +  3) 


x{(2/fc  -  l)(4Ai  +  3)2kG(k)  -  H(k)G(k  -  1)} 


(rt)2fc_2r2 


(|)fc 

k\ 


+E 

k=2 


m 

H(k) 


i  4A1  T  1 

(2A:  +  l)(2fc-l)^  +  F(fc)^±I 
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x2  k 
Using 


(I)* 


i 


A:!  J  ( rt ) 


1 


'2k— 2 


+ 


a 


2H{1) 


1 


A\f\  (  -,0  )  +  A2f3[  — , 0  )  —  Ai  -  A2^. 


-5F(1)  +  ^ 


1 


(3.24)  and  (3.26)  to  obtain  the  following  expres¬ 
sion: 


a 

IT 


-1 


1  +  T 


1  —  T 


G(r) 


a 


2  (r  +  l)  +  l 


we  finally  obtain 


[R(s,t) 

+S(s,t)]cIt  =  1,  —  1  <  s  <  1.  (4.1) 


+  £ 


2{2k  +  l)2 


_2  1+u  1 

7T  ’  3  r3t2  '  4—^  H(k)  I  (A;  +  l)r2 

k= 1 


x  U-ffc  +  l)^)  - 


1 


1 


t2  4(4A;  +  3) 

x2kG{k)  -  H(k)G(k  -  1)}  j  +  A;F(A;)<{  (2A:  +  2) 


The  numerical  solution  technique  is  based  on 
the  collocation  scheme  for  the  solution  of  singu¬ 
lar  integral  equations  given  by  Erdogan,  Gupta, 
and  Cook  [11].  0  This  amounts  to  applying  a 
{(2A;  —  l)(4fc  +  3)  Gaussian  quadrature  formula  for  approximating 
the  integral  of  a  function  /(r)  with  weight  func- 


'{< 


tion  [(l+r)/(l—  t)]2  on  the  interval  [-1,1].  Thus, 
letting  n  be  the  number  of  quadrature  points, 


(2fc-l)!!\2  1 

{2k)\\  J  ( rt)2kr 


i: 


1  +T 
1  —  T 


f(r)dr  = 


2n 


2n  +  l 


5 ~2{1  +  Tk)f(rk ), 

(4.2) 


k= 1 


The  above  equation  completely  agrees  with  the  where 
result  by  Atsumi  and  Shin  do. 

A  quantity  of  physical  interest  is  the  stress  in¬ 
tensity  factor  which  is  given  as 

K=  lim  a/2 (r  -  7 )azz(r,  0). 


Tfc  =  COS 


2k  -  1 
2n  +  1 


7 r, 


k  =  1, . . .  n.  (4.3) 


r— >7“* 

We  choose  p(r)  =  po  and  put  7  —  1  =  0.  We  let 

r  =  |(S  +  l)  +  l,  t=  |(t  +  1)  +  1,  (3.24) 

and  in  order  to  facilitate  numerical  analysis,  as¬ 
sume  g(t )  to  have  the  following  form: 

g(t)  =  Po{t  -  1)5(7  -  t)~*G(t).  (3.25) 

With  the  aid  of  (3.24),  g{r)  can  be  rewritten  as 


The  solution  of  the  integral  equation  is  ob¬ 
tained  by  choosing  the  collocation  points: 

2m 


Si  =  cos 


i  =  1, 


,n, 


2n  +  1 

and  solving  the  matrix  system  for  G *(77.)  : 

2n  +  l 


(4.4) 


k= 1 


where 


2  a 


g{r)  =  p0G(t) 


1  +  r 
1  —  r 


(3.26) 


y>(Sj,7fc)  +  S(Sj,Tfc)]G*(Tfc)  = 

j  =  l,...,n,  (4.5) 

G(rfc)  =  n  1  \  1  11  *  (4-6) 


The  stress  intensity  factors  K  can  therefore  be 
expressed  in  terms  of  G(t)  as 

K/po  =  V^Gfr),  (3.27) 

or  in  terms  of  the  quantity  actually  calculated 

K/P0^l  =  >/2G(7).  (3.28) 


(1  +  Tfe)[f  (Tfe  +  1)  +  1] 

5. Numerical  results  and  consideration 


Numerical  calculations  have  been  carried  out 
for  v  =  0.3.  The  values  of  normalized  stress 
intensity  factor  K/po^/a  versus  a  are  shown  in 
Fig.  1-3  for  various  values  of  5. 

Fig.l  shows  the  variation  of  K/po^/a  with  re¬ 
spect  to  a  when  <5  =  0. 

This  figure  shows  that  as  a  increases,  S.I.F. 
substitutions  are  made  by  the  application  of  decreases  steadily. 


4. Numerical  analysis. 

In  order  to  obtain  numerical  solution  of  (3.23), 
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Fig. 2  and  3  deal  with  the  cases  when  5  =  0.3 
and  5  =  0.5,  respectively.  Here  we  omit  units. 
We  can  see  that  the  trend  is  similar.  Theoreti¬ 
cally  the  infinite  series  involved  converges  when 
5  <  1  by  comparison  test.  However,  because  of 


the  overflow,  computations  could  not  be  accom¬ 
plished  beyond  the  values  6  >  0.5.  And  we  found 
that  the  variation  of  SIF  is  very  small  with  re¬ 
spect  to  the  variation  of  6. 


o 


a 


Fig.  1  Variation  of  s.i.f.  when  delta=0 


a 


Fig, 2.  variation  of  s.i.f.  when  delta=.  1 
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a 


Fig. 3  Variation  of  s.i.f,  when  deita=  5 


6.  Penny-shaped  crack. 

In  this  section  we  are  concerned  with  a  penny¬ 
shaped  crack  in  a  convex  lens  shaped  elastic  ma¬ 
terial.  The  problem  of  determining  the  distri¬ 
bution  of  stress  in  an  elastic  sphere  containing 
a  penny-shaped  crack  or  the  mixed  boundary 
value  problems  concerning  a  spherical  boundary 
has  been  investigated  by  several  researchers.  Sri- 
vastava  and  Dwivedi  [12]  considered  the  prob¬ 
lem  of  a  penny-shaped  crack  in  an  elastic  sphere, 
whereas  Dhaliwal  et  al.  [13]  solved  the  problem  of 
a  penny-shaped  crack  in  a  sphere  embedded  in  an 
infinite  medium.  On  the  other  hand,  Srivastav 
and  Narain  [14]  investigated  the  mixed  boundary 
value  problem  of  torsion  of  a  hemisphere. 

7.  Formulation  of  problem  and  reduction  to 
a  Fredholm  integral  equation  of  the  sec¬ 
ond  kind.  0  We  employ  cylindrical  coordinates 
(r,  0,  z)  with  the  plane  z  =  0  coinciding  the  plane 
of  the  penny  shaped  crack.  The  center  of  the 
crack  is  located  at  (r,z)=(0,0).  As  before,  spher¬ 
ical  coordinates  ( p ,  6 ,  cj))  are  connected  with  the 
cylindrical  coordinates  by 

z  =  pcos6 ,  r  =  psin9. 

Spherical  coordinates  (C,  $,  </) )  whose  origin  is  at 
z  =  —6,  r  =  0,  and  is  the  center  of  the  upper 
spherical  surface  of  the  convex  elastic  body,  is 


also  used.  The  elastic  body  is  symmetrical  with 
respect  to  the  plane  z  =  0. 

The  crack  occupies  the  region  z  =  0,  0  < 
r  <  1.  The  radius  of  the  spherical  portion  is 
Co  =  y/72  +  <52.  The  boundary  conditions  are: 

On  the  plane  z  =  0,  we  want  the  continuity  of 
the  shear  stress,  and  the  normal  displacement: 

uz(r,  0+)  —  uz(r,0~)  =  0,  1  <  r  <  7,  (7.1) 

arz(r,0+)-arz{r,0~)  =  0,  0  <  r  <  7.  (7.2) 

And  the  crack  is  subjected  to  a  known  pressure 
P(r),  be., 

Vzz(r,0+)  =  -p(r),  0  <  r  <  1.  (7.3) 

On  the  surface  of  the  spherical  portion,  stresses 
are  zero: 

o'CcCCo,  V)  =  0,  (7.4) 

*C*(Co,0)  =  O.  (7.5) 

We  can  make  use  of  (2.6)-(2.9)  also,  for  the 
present  case.  The  functions  q and  for  the 
regions  z  >  0  and  z  <  0,  respectively,  are  chosen 
as  follows: 

ro o 

0(1 )(r ,  z)  =  (2u  -  1)  /  C_1^(C)  Jo(COe~^VC 

Jo 
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oo 

+  y^janpnPn(cos  6),  (7.6) 

72—0 

roo 

4>(2){r,  z)  =  (2v  -  1)  /  £-1A(£)  J0(£r)etzd£, 

Jo 


By  substituting  (7.11)  into  (7.12),  it  reduces 

2  00 

g{t) - V^(— 1)n{2nt2n_1a2n 

7T  z ' 

72—0 

+ab2n+1t2n+1}  =  h(t),  0  <  r  <  1,  (7.13) 


-^an(-l)VPn(cos0).  (7.7) 

n= 0 

Here  the  superscripts  (1)  and  (2)  are  taken  for 
the  region  z  >  0  and  z  <  0,  respectively.  The 
functions  ip^  and  ip^  are  chosen  as  follows: 

roo 

ip<'1\r,z)=  A(^)MM)MzdP, 

Jo 

oo 

+  E  bnpnPn(cosO),  (7.8) 

72—0 


roo 

M 2\r,z)  =  -  A(£)J0(£r)etzd£, 

Jo 

oo 

+  ^2bn(-0npnPn(  cos  6).  (7.9) 

n=0 

Then  we  can  immediately  satisfy  condition  (7.2) 
by  these  choice  of  functions  (7.6)-(7.9). 

Now  the  condition  (7.1)  requires 

A(OM*;r)dZ  =  0,  r>  1.  (7.10) 


Equation  (7.10)  is  automatically  satisfied  by 
setting 

MO  =  [  g(t)sm(£t)dt,  g(0)=0.  (7.11) 

Jo 

Then  from  the  boundary  condition  (7.3),  we  ob¬ 
tain 


M(0 MM)<%  -  Y,  a2n(2n)2P2n(0)r2n~2 

n=0 


“2  Wl^L+lCO)^2"  =  ~p(r), 

n= 0 

0<r<l,  (7.12) 

where  the  prime  indicates  the  differentiation 
with  respect  to  the  argument.  0 


where 


The  solution  will  be  complete,  if  the  condi¬ 
tions  on  the  surface  of  the  spherical  portion  are 
satisfied. 

8. Conditions  on  the  surface  of  the  sphere. 

Equation  (7.13)  gives  one  relation  connecting 
unknown  coefficients  an  and  bn.  The  stress  com¬ 
ponents  besides  (2.8)  and  (2.9)  which  are  needed 
for  the  present  analysis  are  given  by  (3.1)  and 
(3.2).  To  satisfy  boundary  conditions  on  the 
spherical  surface,  it  is  needed  to  represent  <p,  ip 
in  (7.6)-(7.9)  in  terms  of  £,  d  variables.  To  do  so 
we  utilize  the  following  formula  in  the  Appendix 
2.  An  expression  useful  for  the  present  analysis 
is  the  following 


Pn(cos6)pn  =  E  (  I  )  Pk(cosd)Ck(-S)n~k 


Thus 


(8.1) 


OO  OO 


E  anPn( COS  0)pn  =  E  a?l  E  (  t  )  Pk(COS  &) 


n= 0 


n= 0  k= 0 


n\ 


x(k(-5)n~k  =  £pj(co80)C jAj,  (8.2) 

3=0 

where 

Also 

r)  r°° 

~¥tJo  C1mMZr)e-tzdS 

r  1  roo 

=  ~  g(t)  /  MM)  cos (£t)e~Zzd£dt.  (8.4) 
Jo  Jo 


^  .  ^.an(-On~J-  (8-3) 
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The  inner  integral  of  the  above  equation  is 

POO 

/  Jo(£r)cos(£t)e_^d£ 


The  inner  integral  on  the  right-hand  side  of 
(8.10)  is 


/  o 


=  3F l 


POO 

/  MCr)e~^z+itUc 

Jo 


1 


=  &- 


y/r2  +  (z  +  it)2 

1 

yjr2  +  z2  +  2  zit  —  t2 

1 


E 

n= 0 


/  MJr)e-^z+lt)d^ 

Jo 

P2n+ l(cosfl)/  j\n^2n+l 

^271/~\~  2 


(8.11) 


To  express  (8.11)  in  the  (£,  '&)  spherical  co¬ 
ordinates,  we  use  the  following  formula  in  the 
Appendix  1. 


TWi(cos  9) 


p^l-2cos9(=f)  +  {=^)2 

1  00  /  t  \  2n 

El 


p 


2n+l 


E 


(8.5) 


x5 


k= 0 

kP2n+l+k{cOSd) 


(2n  +  l  +  k)\ 
k\(2n  +  1)! 


(8.12) 


^  2 77.  2 /c 

If  we  use  (8.8)  and  (8.12),  we  get  i/a1)  in  spherical 
Thus  finally,  from  (8.2),  and  using  the  for-  coordinates  (C,d)  as 
rnula  in  Appendix,  cj)A>  can  be  written  in  terms  oo  ,  ^  >. 

of  spherical  coordinates  (£,  •&)  as  =  E  Pk(cosd)  (  +  Bk(k  )  +  B0, 

k= 0  ^  ^  ' 

(8.6) 


0(1)  =  ^Tfc(cosd) 
where 


k=0 


Afe+ 


i+Afcc* 


(8.13) 


where 


[fc/2] 

=  -(a  -  1)  Y 


k\ 


[(fc-l)/2] 

E 


k\ 


rk—2n—l 


rk— 2n 


n= 0 


(2ra  +  l)!(fc  —  2n  —  1)! 


71=0 


(2n)\(k-2n)\ 


x(-i r 

v  1  Jo  2n  +  l 


x(-l)n  /  5(t)f2n+1df. 


(8.14) 


(8.7)  Then  if  we  substitute  these  values  of  and 
^l1)  given  by  (8.6)  and  (8.13)  into  (3.2),  and 
It  is  also  necessary  to  express  in  terms  of  simplify  the  results  by  using  the  properties  of 
spherical  coordinates  ((,$)■  Now,  as  in  (8.2)  Legendre  polynomials,  from  the  condition  that 

on  the  spherical  surface  C  =  Co5  =  0,  we 
^  oo  obtain  the  following  equation 


Y  bnPn( cos  0)pn  =  Y  Pj (cos  1 (8.8) 

n=0  j= 0 

where 

OO  , 

n\ 


Ak+lk Ct1  -  Pk+2( 0 


k+i  a  —  {k  +  2y 


2k  +  5 


_  k_1k{2a  -  k  -  1)  fc  +  3 

~B^o  notl1  -  Tfe+4 


bn(-5)n~P 


n=j 


,  .  2A:  + 1  Cn 

(8.9) 

k+2  (fc  +  3)(fc  +  2  +  2a)  a  —  (A:  +  l)2 


0  We  first  express  following  integral  in  (p,  9)  co¬ 
ordinates 

POO  P 1 

/  -4(6 Jo(COe_5z^  =  /  s(*)di 

Jo  Jo 

POO 

x  /  sin(£t) Jo(^r)e~^zd^.  (8.10) 

Jo 


/-k+4 

’0 


2/c  +  5 


=  0. 


^+2  2k  +  1 

(8.15) 


Likewise,  from  the  condition  cr^(Co;d)  =  0, 
we  get  the  following  equation, 


—Ak+ik(k  +  1)C0 


fc-i 
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+-B/C+2C0 

+-®fcCo 


fc+i  (&  +  2){2  —  a  —  (k  —  1  )(k  +  2)} 
2k +  5 

fc_i  k(k  +  1)(2  a  —  k  —  1) 


(2k  +  1) 
(k  +  2)  (A:  +  3) 


I  (7;  I  i)(k  I  5)  G(fc)  I  G(fc  +  2) 
+(A  +  4j(/e  +  5)H(fc  +  2)  +  2jfc  +  5  j, 

G(fc) 


JV(ife)  = 


(2k  +  7)  -  1 


c, 


k+ 4 


-T 


(fc  +  2)(k  +  3)(A;  +  2  +  2a) 
k+2  Co+4(2  k  +  5) 

(A:  +  l){(/c  +  4)(fc  +  1)  +  a  -  2} 


H(k  +  2) 

(k  +  5)I(k  +  2) 

(2k  +  9)(k  +  2)(k  +  3)J 


with 


C0fc+2(2fc+l) 

=  0. 

Thus  if  we  multiply  (8.15)  by  k  +  1  and  add  Pen^ix  we  have 
(8.16),  from  the  resulting  equation  we  find 


H(k)  =  -a(2k  +  3)  +  2  (k  +  2  )(k  +  3), 

G(k)  =  2—a—(k+l)(k+4:),  I(k)  =  —a+(k+ 3)2. 
(8.16)  Therefore  using  the  formula  in  Theorem  B  of  Ap- 


Bk+2  k+1 


(2k +  5) 


Co  {— ot(2k  +  3)  +  2(k  +  2)(k  +  3)} 

(k  +  3)(2k  +  3) 


s-k+A 

SO 


‘f’fc+i 


+^fc+2 


(2k  +  3)(k  +  3)(k  +  2  +  2a) 


C0fc+4(2fc  +  5) 


vp , 

+  77+2  (^’  +  1)(fc  +  3) 
So 


Using  the  relation 


Bk+ 2  is 


Bk+2  k+1 


,  a  -  1 

k+1  ~  ~  k  +  2^fc+2’ 

(2fc  +  5) 


Cq  {— a(2k  +  3)  +  2  (k  +  2)(fc  +  3)} 

(2k  +  3)(k  +  3)((k  +  3)2  -  a) 


^k+2 


Co+\2k  +  5)(k  +  2) 


+  ~7kT' 2 +  1)(^c  +  3) 

so 


From  (8.15)  we  have 


A.+3  =  +  SgM(t)  +  5±ijv(t). 


so 


so 


where 


M(fc)  = 


L(k )  =  (k  +  l)(k  +  3) 

1 


So 

(2a- k-  3) 


k  +  2 


(k  +  3) 


W(fc)  ’ 

(2k  +  3)  (2a  —  k  —  3)I(k) 
(2k  +  5  )H(k) 


A  =  ^(-l)n2nt2re-1a2n 


n= 1 


^(-l)n2nt2n-1 


k\Ak5k 


8‘ln(2n)\  ~^n  (k  —  2n)\ 


E 


k\(-l)r 


A2nl 


(8.17) 


n=  1 

OO  r[&/2] 

^i(2n-l)\(k-2n)\\5)  t 

OO  OO 

xAk5k  =  X  fk(t)Ak+ 3  =  X  /fc(i) 


fc=-i 


fc=-i 


**  m  +  StlM(t)  + 


/-2k+3 

LSo 


/■ 2fc+5 ' 

SO 


/-2fe+7 
SO 


(8.19) 


In  the  above  equation  ^>k  =  0,  if  k  <  0  and 
[(fc+3)/2] 

fk(t)  =  Sk+3  J] 


(fc  +  3)!(-l)' 


n= 1 


(2n-  l)!(fc  +  3-2ra)! 


t\2nl 


5/  t 


(8.18)  If  we  substitute  the  values  of  \Pfc  into  (8. 19), it 
reduces  to 

A  =  I  g(u)Qi(t,u)du, 


where 


Qi(t,u)  =  X]  fk(t) 


+ 


hk+ 2(u) 

Ci 


2fc+5 

0 


fc=-l 


Af(fc)  + 


h-Mm 


LC( 


hk+A(u) 
/-2k+7 
S( 


JV(fe) 


0 
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In  the  above  equation  hk{u )  =  0,  if  k  <  0  and 
[(*- 1)/2] 

hk(u)  =  6k  J2 


(-1  )nfc! 


=  C2i  1  -  2 cos i9-  + 


t'o  (2n+l)!(fc-2n-l)! 


1 


*'$) 


2n+l 


<5  i  ft\2 
Pn( cosd)  (-l)n  dn  ( 1 


p  CA/l  -  2cosi?|  +  (|)~  71=0 


X>(costf) 


£n+l  ’ 


Also,  using  Theorem  B  of  Appendix  4,  we  get 

^.2n+l 


0n+l 


n! 


^(-l)"62r!.+1t2"+1  =  £(-!)"- 


n= 0 


n=0 

00  tit -»  cfc 


(2n  +  l)!<52n+1 


E 

/c=0 

oo 


(  l)n  /P*(cos0)\.fc 


n!  Szn  V  Cfc+1 


E 


jfe!sfc<y* 


r  fc— l 
OO  |-L  2 


fc=2n+l 

(— l)nfc! 


(fc  —  2n  —  1)! 


(_l)n  gn+fc  /  j\  (-l)fc 

^  n!  0^+fcVC/  fc! 


/  ,  \  2n+l' 

\  \  _ l-1!  _ (  -) 

t^o  (2n  +  l)!(fc  —  2n  —  1)!  \SJ 


Bu5k 


k= 0 
oo 

E 

k= 0 


n  +  fc  \  Pn+k( cos  i?)  fc 
A:  y  £n+fc+i 


F,  hk+2(t)Bk+2 

k=- 1 

2/ +  5 


Appendix  2.  Proof  of  (3.16). 


>/c+2 


fe=_X  C+iH(fc)  LCo 

^fc+2  (fc  +  3)  (2 k  +  3) 


(k  +  l)(k  +  3)  2-k  pn  Pn{cos  0)  =  /  (z  +  ix  cos  u  +  iy  sin  u)ndu 


+ 


0}~+4  2k  +  5 


/(fc) 


=  /  (— <5  +  Z  +  ix  cos  u  +  iy  sin  u)ndu 


-i 


=  /  g(u)Vt2(t,u)du, 

Jo 


E 

fc=0 


n! 


(n  —  k)\k\ 


M) 


n—k 


0  where 


„ .  ,  .  I  /  X  2fc  +  5  |"/lfe('u) 

n2(t,u)=  Y,  hk+2(t)  -Mik  +  i) 

k=~ i  So  rtW  L  So 

x  (fc  +  3)  +  MM7(A) 


x  /  (Z  +  ixcosu  +  iy  sin  u)kdu 

J  —TV 

n  . 

n! 


So 


2k  +  5 


2»E 


(-r'c'aw. 


fc=0 


(n  —  fc)!/! 


In  the  above  equation  /ifc(«)  =  0,  if  fc  <  0.  Thus  Appendix  3.  Theorem  A.  If  the  equation 


(7.13)  reduces  to  the  following  Fredholm  integral 
equation  of  the  second  kind 

9(t)+  [  g{u)K(t,u)du  =  h(t), 


a  =  E 


fc! 


^  (fc  —  *)!*! 

7=0  V  ' 


is  solved  for  aPs,  it  will  be  written  as 


where 


K(t,  u )  = - {f2i(t,  u )  +  o&l2(t,  it)}. 

7T 


Appendix  1.  Proof  of  (3.3).  Since 

p2  =  (C  cos  i?  —  5)2  +  (C  sin  i9)2 


di  — 


^  a 


k= 0 


(i  —  fc)!fc! 


5fc(-l) 


i—k 


Proof.  We  prove  it  by  mathematical  induction. 
Suppose  it  is  true  for  i.  we  will  show  that  it  is 
also  true  for  i  +  1.  Suppose  //+i  is  given  by 

d  —  (*  +  !)!  _ 

-£*2+1  —  /  j  /  •  |  i  7  \  |  /  j  ®k  —  1 


k= 0 


(*  +  1  —  fc)!fc! 
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(*  +  !)• 

+  77TT,ak- 


Proof.  Let 


k= o 


(*  +  l-ife)!fc! 


Thus 


ai+i  -  Bi+ 1  ^2  /  ■ 


(*  +  !)! 


k= 0 

k\ 


3= 0 


3  Kk~j) 


(i  +  1  -  Jfc)!fc! 


then 

and 


J2an(~S)nXn  =  /(*), 

n=0 

an(-«5)"n!  =  /W(0), 


-  5i+i  2^  Bi  41  2^  a 


(-l)fc-i 


J=0 


J! 


k=j 


(i  +  1  —  k)\(k  —  j)\ ' 


n=k 

From  Taylor’s  series 


(Al) 


The  inner  summation  in  the  above  equation  can 
written  as,  by  changing  the  variable  k  —  j  =  m  Thus 


fc=o 


k\ 


E?7 


_1  i k-j 


(-1) 


k=j 

i~3 


(i  +  1  -  k)\(k-  j)\ 

(-l)m 


>M)nn!  =  /(")( 0)  =  £ 


k= 0 


/(fc)( i)  dn 
k\  dxn 


(i-d‘ 


x=0 


Et 7 


V  /W(l)  r  i  \k—n  =  V  ^(-1)fc_n(-‘J)fcfc! 


— J  (i  +  1  —  m  —  j)\m\ 

m= 0  v  J  ' 


k=n 

Therefore 


k=n 


(k  —  n)\ 


(' i~j  +  1)! ' 


CLn  — 


(-l)m(i-j  +  l)! 


since 
i-j+l 

y  x 

(i  +  1  —  m  —  j)\m\ 
m= o  v  JJ 

Then  (A.l)  is  equal  to 


E 

k=n 


k\ 


n\(k  —  n) 


Bk5 


k—n 


□ 


=  (i-i) 


i-j+ 1  _ 


=  0. 


CH+ 1  — 


Bi+i  +  'y 


U_1  v+i-i 


i=o 


Bj(i  +  1)!(  1) 

j!(i  +  1  -  j)! 


i+1 

E 

i=o 


u_i  ++i-i 


Bj(i  +  !)!(-!) 

j!(i  +  1  -  j)! 


□ 


Appendix  4.  Theorem  B.  If  the  equation 


b*  =  E 


n=k 


n 


(-S)n~kan, 


is  solved  for  a^s,  it  will  be  written  as 

OO 

E 


dr).  — 


k=n 


n  ■ 
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